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>^' Abstract 

c^ : 

So far there is no closed formula for relative entropy of entanglement of arbitrary two-qubit 

^O I states. In this paper we present a method, which guarantees the derivation of the relative en- 

tropy of entanglement for most states that have z-directional Bloch vectors. It is shown that the 

Q_i' closest separable states for those states also have 2:-directional Bloch vectors though there are few 

rj . exceptions. 
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Research into entanglement of quantum states has long history from the very beginning 
of quantum mechanics l|, |2|. At that time the main motivation for the study of entangle- 
ment was to explore the non-local property of quantum mechanics. Still, this issue is not 
completely understood. Recent study on the entanglement is mainly due to its role as a 
physical resource for the various quantum information processing such as teleportation|3|, 
quantum cryptography [J], and speed-up of quantum computer |5|]. 

In order to quantify how much a given quantum state is entangled many entanglement 
measures were invented for last two decades. Among them the most important measure 



u 



seems to be the distillable entanglement [6[] because it measures how a given quantum state is 
useful in the real quantum information processing with overcoming the effect of the noises via 
the purification protocol. In spite of its importance the analytically derivational technique 
for it even in the relatively simple quantum system is not known. In fact, in order to compute 
the distillable entanglement we should find an optimal purification protocol. However, it 
is nontrivial problem to find the optimal protocol except very rare cases. In this reason 
many people tried to find more analytically tractable entanglement measures which may be 
able to provide an information on the tight upper bound of the distillable entanglement. 
The representatives constructed in this reason are entanglement of formation (EOF) 161 and 

nn 

relative entropy of entanglement (REE)[7|, |8[. 

About a decade ago Woottersfo!] found how to compute the EOF for arbitrary two-qubit 
states. Although still we do not have closed formula of EOF for higher- dimensional quantum 
system, the Wootters' result has great impact in the study of entanglement. One of the 
example for an application of the Wootters' result is to examine the role of the quantum 
entanglement in the complex quantum system such as bio-system [10]. Another direction of 
application is to use the Wootters' result to find a truly multipartite entanglement measure. 
In this way, the three-tangle, measure for the genuine tripartite entanglement, was invented 



inRef.(n|. 



On the contrary, still we do not have closed formula of the REE even for the two-qubit 



states[l2|]. In order to understand the distillable entanglement more profoundly, therefore, 
it is worthwhile to investigate the properties of the REE for the various two-qubit states. In 
this paper we would like to examine the REE for the states, which have ^-directional Bloch 
vectors. We present three theorems in the following, which guarantees that the REE for the 
most such states can be computed analytically or, at least, numerically. 



The REE for state p is defined as 



Er{p) = mill S{p\ \a) = min ti [p\n p — p In a] 



(1) 



where D is a set of positive partial transpose (PPT) states. For various properties of the 
16| . If our concern is restricted into the two-qubit state, it is possible to 



REE see Ref. 
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regard P as a set of the separable states, because there is no bound entangled state in the 
two-qubit Hilbert space. The separable state a in Eq. ([1]) is called the closest separable state 
(CSS) of p. In order for the separable state a to be CSS of some entangled states it should 
be edge state in the set T>, which means that the smallest eigenvalue of a^ is zero, where 
the superscript F denotes partial transposition^. 

Although the definition of the REE is comparatively simple, the analytic computation 
of it is highly difficult problem even for the most simple two-qubit case (see chapter 8 
of Ref.[l2|). Since the REE can be straightforwardly computed provided that the CSS is 
derived, this means that finding a CSS of the given entangled state is very difficult. Recently, 
however, the authors in Ref. 17|] analyzed the converse procedure. When the edge separable 
state TT is full-rank, they have presented a method for deriving the entangled state p, whose 
CSS is TT. Still, however, finding a CSS for the arbitrary entangled state p is an unsolved 
problem. 



In order to explore the issue for finding CSS or REE authors in Ref. 18| approached the 
problem from the geometrical point of view. To explain the main results of Ref. [l8| briefiy 
it is convenient to express the given entangled state p in a form 
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P = -A 



r ■ cr 



I + I' 



S ■ (T 






CTi 



cr^ 



(2) 



where cr is usual Pauli matrices. The vectors r and s are Bloch vectors for each qubit and 
the tensor Qij represents a correlation between qubits. Since appropriate local-unitary (LU) 
transformation for each qubit can make the correlation tensor Qij to be diagonal, it is more 
convenient, without loss of generality, to express p as 
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I®I + r-(7®I + I®s-cr + ^g^ 



P = l 



i\^ n 



a„ 



n=l 



(3) 



^ The converse of this statement, i.e. if a is an edge state in D, there exist entangled states whose CSS are 
a, is not generally true 



For example, for the four Bell states 
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|/3i 



V2 



(|00) + |11)) 



1/32) 



v^ 



(|00)-|11)) 



(4) 



\(3,) = —{\0l) + \10)) 



|/34) = -^(|01)-|10)) 



the Bloch vectors r and s are vanishing and the corresponding correlation vectors become 

91 = (1,-1,1) 92 = (-1,1,1) 93 = (1,1,-1) 94 = (-1,-1,-1). (5) 

In Ref. 



18j it was shown that if p is one of Bell-diagonal, generalized Vedral-Plenio (VP) 



and generalized Horodecki states, its CSS is 

/(g)/ + r-cr(g)/ + / 



TC 



S ■ (J 



+ J2^r. 



n=l 



CTr, 



CTr. 



(6) 



The correlation vector of tt, 7, can be computed from a fact that the straight line in the 
correlation vector space, which connects 7 = (7x,7j/,7^) and 9 = {gx,gy,gz) passes through 
one of Eq.(|5]), which is the nearest one from 9. Since this fact with the edge state criterion 
uniquely determines the correlation vector 7 of the CSS, it is straightforward to compute the 
REE for the Bell-diagonal, generalized VP and generalized Horodecki states. For example, 
let us choose the Bell-diagonal, VP and Horodecki states as following: 

Pb = Ai|/3i)(/3i|+A2|/32)(/32|+A3|/33)(/33|+A4|/34)(/34| (7) 

(max(Ai,A2,A3,A4) = A3) 
Pvp = Ai|/33)(/33|+A2|01)(01|+A3|10)(10| 
Ph = Ai|/33)(/33|+A2|00)(00|+A3|ll)(ll|. 



Following Ref. 18| it is easy to show that the corresponding CSS for these states are 
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(Ai + 2A2)(Ai + 2A3) 



+ A2 |01)(01|+hr + ^3 |10)(10| 
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and their REE become 



|/33)(/33| + ^H^|00)(00|+i^ 



2A3)^ 



11)(11| 
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-^ryPvp) 
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In 2 
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H[^ + X,]~H{A) 



A=l 
2 



Er{p 



H 



Ai In Ai + A2 In A2 + A3 In A3 + 2if ( — + A2 



Ai 



1 + \/a? + (A2 - A3)^ 
- Ailn2 



where H{p) = —pimp — (1 — p) ln(l — p). It is worthwhile noting that Er{pvp) and Er^pn) 
are invariant under the exchange of A2 and A3. In fact, one can conjecture this symmetry 
from the physical point of view. 

In this paper we would like to examine the REE for the two qubit states, whose Bloch 
vectors r and s are ^-directional. Thus, we assume r = (0,0, r) and s = (0,0, s). For more 
simplicity we assume that the first two components of the correlation vector g are identical, 
i.e. gx = Qy. Then, the quantum state p can be written as 

f Ai \ 

A2 De'^ 

De-'^ A3 

\ A4/ 



P 



(10) 



where 



Ai 



l + r + s + Qz 



A, 



1 + r - s- g^ 



(11) 
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We also impose 



1-r + s - g^ 



A, 



r - s + g^ 



D 



9x 



2 cos 99 



> 0. 



D^ > AiAi 



(12) 



TT 



(13) 



to require that p is an entangled state. 

Now we conjecture that the CSS of p is of a form 

/ri 0\ 

ra ye'^ 

ye-'^P n 

\ r4 / 

with y = y/rir^ < ^/r^^Xz- In the following we will show that most entangled states of the 
form flTOj) have really their CSS as the form flT3l) . However, for extremely asymmetric states 
we will show that our conjecture is not true. 



If TT is really the CSS of p, the following coupled equations should be satisfied [17|: 
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Ta 
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ri +r4 
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ri +r4 
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(14a) 
(14b) 



27^1 7^/1 

^2 + x-( \ VT-o [2'^i'^4^ + (r2 - r3)(r2£ - z)] = A^ (14c) 

27*1 7^/1 

^3 + a; 7 ^ r^ [2rir4£ - (ra - rg) (rgf - 2;)] = A3 (14d) 

7/ 

y + a; 7 — 7 — x-^ [2rir4(r2 + rs)^ + (r2 - rs)^^;] = D, (14e) 



where a: is a positive parameter and 



'^2 + ^^3 + 2; 



^ 2 



-2 = \/(r2 — r-i)^ + 4rir4 £ = In . (15) 

r2^r-i- z 

In this case one can show after tedious calculation that the REE of p becomes 

Er{p) = tr(p In p) — tr(p In tt) (16) 

= {Ai In Ai + Ai In A^ + ^4+ In A+ + A_ In A_) 

(.. _^ . 1 , A2 + A3 , , (A2-A3)(r2-r3)+4Dy 
- I y4i In ri + y44 In r^ H In(r2r3 - rir^) H TT^i 

where 

[{A-2 + As) ± ^(^2-^3)2 + 4^2] . (17) 

Now we present the following three theorems, which provide the REE and CSS of the 
entangled state p given in Eq. lfTOl) . 

Theorem 1. If Ai = ^4 = 0, Er{p) becomes 

Er{p)=H{A2)-H{A+). 

Proof. If Ai = A4 = 0, Eq.( 114ap and Eq.( 114b[) give solutions ri = r^ = e, where e is an 
infinitesimal positive parameter, which will be taken to be zero after calculation. Then, the 
remaining equations in Eq. lTT^ eventually generate the following solutions: 

\A2-A3\ mm(A2,v43) 
Therefore, CSS vr in this case is 

7r = A2|01)(01|+A3|10)(10|. (19) 

Making use of Eq. (ll6p it is straightforward to compute the REE, which completes the proof. 
As an example of theorem 1 let us consider 

p = p\4^)m+qi\01){01\+q2m{l0\ (20) 



where p + gi + gs = 1 and |^) = a|01) + /3|10) (lap + |/3|2 = 1). Then the CSS of p is 

TT = (p|a|2 + gi)|01)(01| + (p|/3|2 + g2)|10)(10| (21) 

and the corresponding REE is 

Er{p) = H {p\a\^ + gi) - H{A+) (22) 

where 



2 L 



1 ± ^f^^+I^^^^HM^¥^\m + l^^^^\ ■ (23) 

When a = l3 = 1/v^, it is easy to show that Eq.(!22l) reduces to the second equation of 



Eq.(l9]) when Ai = p, X2 = Qi and A3 = q^. 

Theorem 2. If both Ai and A4 are not zero, and A2 = A^, the REE of p becomes 



Er{p) = n,-n2 (24) 



where 



In Eq.(l25 



^1 = Ai In Ai + Ai In A^ + {A2 + £>) ln(v42 + D) + {A2 - £>) ln(A2 - D) (25) 

r22 = Ai In ri + ^4 In r4 + y42 In(r2 — titA + -D In . 

r2-y 



n = j [2 Ai ( Ai + A2) ( Ai + A2 + Ai) - D^ {A, -Ai) + A] (26) 

r4 = 4; [2^4(^2 + Ai){Ai + A2 + A4) + D\Ai - A^) + A] 

r 

r2 = - [2(Ai + A2)(A2 + A4)(Ai + A2 + A4) - D\Ai + 2^2 + A4) - A] 
r 

where y = a/t^T^ and 

F = 2{Ai + A2 + A4 + D){Ai + A2 + Ai-D) (27) 

A = D^/D^{A^ - A^y + AAiA^iA^ + A2)iA2 + A^). 

Remark: Under Ai -v^ A4, r2 is invariant and, ri and r4 are changed into each other. 
This fact indicates that E^{p) is invariant under Ai ^ A4. The appearance of this symmetry 
is plausible from the physical point of view. 



Proof. Since both Ai and A4 are not zero, Eq. fll4ap and Eq.f ll4bp enable us to express 
r^ and x in terms of ri as follows: 

n = n- {A, -A,) x= ^^^- ^^)^^^ + ^^) _ (28) 

rir4 



Since ^2 = ^43, Eq.f ll4cp and Eq. (ll4dp imply r2 = r^. Then inserting r2 = r^ and Eq. fl28|) 
into Eq. (ll4cp . one can express r2 in terms of ri as follows: 

r2 = -ri + {Ai + A2). (29) 

In fact, Eq. (l29|) can be derived from a normalization ri + 2r2 + r4 = 1. Finally, we consider 
Eq.( ll4el) . which reduces to 

y2-Dy + r2(ri-Ai) = 0. (30) 

Thus, one can express y in terms of ri as a form 

1 



D ± ^ D^ - Ar^iri - Ar) . (31) 

z, L 

Since y^ = rir^, one can compute ri from Eq. (!3T]) . which is 



y=2 



n = ^ [2A,{A, + A2){A, + A2 + A4) - D\A, - A^) ± A] . (32) 

Therefore, one can easily compute r2 and r^ by making use of Eq. (!28|) and Eq. (!29|) . The 
undetermined sign can be fixed by Eq.( l30l) . Then, Eq.( TT6l) completes a proof of theorem 2. 
As an example of theorem 2 let us consider 

P = Pi|/33)(/33|+P2|/34)(/34|+gi|00)(00|+g2|ll)(ll| (33) 

with Pi +P2 + 91 + Q'2 = 1- Then, it is straightforward to show 

2qi{pi+P2 + 2qi){p,+p2 + 2q, + 2^2) - {pi -P2f{qi - 52) + 4A 

ri = (34) 

8(pi + gi + q2)[P2 + gi + 92) 

{pi +P2 + 2gi)(pi +P2 + 2q2){pi +P2 + 2gi + 2^2) - {pi - P2f - 4A 
r2 = 



8(pi + qi + q2){P2 + qi + ^2) 



where 



A= ^ ^ V4gig2(pi +P2 + 2gi)(pi+p2 + 2g2) + (Pi -P2)^(gi -g2)^ (35) 

and r^ is obtained from ri by exchanging gi and g2- Then it is easy to compute the REE 
of p by making use of theorem 2. When p2 = 0, it is also straightforward to show that the 
REE of p reduces to third equation of Eq.dH]) if one identifies Ai = pi, A2 = qi and A3 = ^2- 
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Theorem 3. For other cases the CSS of p can be obtained by solving an equation 



r2 + rj, + z 
r2 + r3- z 



exp 



z{ri - Ai){r2-r3f 



y{D - y)z'^ - 2rir\{r-i - Ai){r2 + r^) 



(36) 



where 



r4 = ri - [Ai - Ai 
1 
4 



(37) 



^3 



(4Ai + 3^2 + A3) - 4ri + Vf 
(4Ai + A2 + 3A3) - 4ri - yf 



and 



r = 16Dy^^ - 8(2^1 + A2 + A^ + 2A^)n + [{A2 - A^f + 8Ai{2A, + A2 + A3)] . (38) 

Remark 1: If Eq.( l37|) and Eq. (!38|) are used, one can make the Ihs and rhs of Eq. (!36|l in 
terms of ri only. Thus, Eq.( l36|) is an equation with only one variable, which can be solved 
analytically or numerically. 

Remark 2: If Eq. fl571) does not provide a solution for some entangled state p, this fact 
indicates that the CSS of p is not of the form flT3|) . In this case, therefore, CSS of p seems 
to have different structure from p. 

Proof. From Eq. (ll4al) and Eq.( ]14b[) one can express r^ and x in terms of ri, which is 
exactly the same with Eq. fl28|) . The remaining equations in Eq. flT^ reduce to 

2z{n - A^){r2 - r^) = i [{r2 - ^2)^' + 2{n - A^) {r2(r2 - r^) + 2^^}] (39a) 

2z{n - Ai){r2 -r3)=i [(A3 - r3)z^ + 2(ri - Ai) {r3(r2 - r^) - 2nn}] (39b) 

2y{D - y)z^ - Arir^^ri - Ai){r2 + rg) 



2z{ri - Ai){r2 - r^] 



(39c) 



r2 -rs 

Since the Ihs of Eq. (l39l) are all identical, the rhs of them should be equal. By equalizing the 
rhs of Eq. (I39ap with the rhs of Eq. (I39cl) one can derive 



(r2 - r3)(v42 - r2) + 2y(D - y) - 2r2{n - Ai) = 0. 



(40) 



Similarly, one can derive 



(r2 - r3)(A3 - r3) - 2y{D - y) + 2r3(ri - Ai 



(41) 



from Eq.f l39bp and Eq. fl39cp . Adding Eq. (HO|) and Eq. (HT]) one can express r2 + r^ in terms 
of ri in a form 

^2 + '^s = 1 - '^i - ^4- (42) 

In fact, Eq. fl42p is a normalization for the CSS vr. Combining Eq. fHT]) and Eq.f l42p one can 
make the following second degree equation 



2rl + [An - (4Ai + A2 + 3A3)] rg + [A3 {(2Ai + A2 + A3) -2n}- 2y{D - y)] = 0, (43) 



which has roots 



1 r 



r-i 



(4Ai + A2 + 3A3) - 4ri ± Vf . 
^ L 

Inserting Eq. f H4|) into Eq. (142 p one can express r2 in terms of ri as a form 

1 



'^ = 1 



(4Ai + 3^2 + A3) - 4ri T Vf 



(44) 



(45) 



The undetermined sign in Eq. fl44p and Eq. fHSj) can be fixed by Eq. fHOl) . Finally, the param- 
eter Ti is determined by Eq.f l39cp . which reduces to Eq. (!36l) . This completes a proof. 

As an example of theorem 3 let us re-consider the model which was considered by Rains 
in Ref.[l9|, where the entangled state is 



'h 











o\ 



Q 45907 7^ 1201 , 49g n. 

^ 90000 150 3750 ~'~ 3600 

n 1201 _|_ 49g 29093 , 7g n, 

^ 3750 "•" 3600 90000 "*" 150 ^ 

Vo i/ 



(46) 



with ^ = l/ln(||). Then Eq.( l36|) directly gives ri = 1/6 and the resulting CSS of p is 



/^ o\ 

6 



TT 



^144 6 

Vo ly 



This is in agreement with Rains' result. 

As a second example of theorem 3 let us consider 



(47) 



P = p|/33)(/33|+gi|01)(01|+g2|10)(10|+g3|00)(00Kg4|ll)(ll| 



(48) 
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with p = 0.66, gi = 0.16, g2 = 0.03, q^ = 0.06 and q^ = 0.09. Then, Eq. (l36|) cannot be 
solved analytically. The numerical calculation shows that the CSS is 

vr = p'|/33)(/33|+g;|01)(01|+g^|10)(10|+g^|00)|(00|+g^|ll)(ll| (49) 

where p' = 0.306933, q[ = 0.252429, q'^ = 0.132241, q'^ = 0.139198 and q'^ = 0.169198. 

Numerical calculation shows that most entangled states of the form (fTOj) have their CSS 
as a form of ( TT3|) . However, there are states whose CSS are not of the form ( IT3l) . For 
example, the state (HSj) with p = 0.66, gi = 0.05, g2 = 0.07, q^ = 0.04 and q^ = 0.18 does 
not have CSS of the form ( 1T3|) . It seems to be interesting to derive a criterion that clarifies 
which entangled states p do not have CSS of the form f fT3|) . 

We have assumed ab initio that the Bloch vectors of p are 2;-directional. In addition, 
we have assumed that the first two components of the correlation vector are equal. These 
assumption are chosen only for simplicity. In the near future we would like to re-visit the 
REE problem for two-qubit states with removing these assumptions as much as possible. 
This may shed light on the explicit derivation for the closed formula of REE in the two-qubit 
system. 
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